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We investigate chiral and deconfinement transitions in the framework of the strong coupling lattice QCD for 
color SU(3) with one species of unrooted staggered fermion at finite temperature and quark chemical potential. 
We take account of the leading order Polyakov loop terms as well as the next-to-next-to-leading order (1/g 4 ) 
fermionic terms of the strong coupling expansion in the effective action. We investigate the Polyakov loop 
effects by comparing two approximation schemes, a Haar measure method (no fluctuation from the mean field) 
and a Weiss mean-field method (with fluctuations). The effective potential is obtained in both cases, and we 
analytically clarify the Polyakov loop contributions to the effective potential. The Polyakov loop is found to 
suppress the chiral condensate and to reduce the chiral transition temperature at /i = 0, and the chiral transition 
temperature roughly reproduces the Monte Carlo results at j3 — 2N c /g 2 < 4. The deconfinement transition is 
found to be the crossover and first order for light {anno < 4 at /3 — 4) and heavy quark masses, respectively. 
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I. INTRODUCTION 

The phase transition in Quantum Chromodynamics (QCD) 
at finite temperature (T) and/or quark chemical potential (/i) 
is attracting much attention in recent years. When thermo- 
dynamic parameters such as T reach to a typical QCD scale, 
abundantly formed hadrons start to overlap with each other 
and matter with color charges would appear, thus a quark- 
gluon plasma (QGP) 01 is expected to emerge at such high 
T and/or fb. The first principle studies based on the lattice 
QCD Monte-Carlo (MC) simulations actually indicate a phase 
transition from hadron phase to QGP around T ~ 160-190 
(MeV) 131. Various experimental observations and theoreti- 
cal arguments imply the formation of a strongly coupled QGP 
in heavy-ion collisions at RHIC. This extreme form of matter 
will be extensively investigated in LHC-ALICE experiments 
at CERN, where we can simulate the matter evolution in the 
early universe. We also expect the transition to quark mat- 
ter when the baryon chemical potential N c /j, overcomes the 
nucleon mass energy N c fi > Mm- The color superconduc- 
tor could appear at high density due to the attraction between 
quarks in the color antisymmetric channel JH, and confined 
dense quark (quarkyonic) matter is also expected from large 
N c arguments |4|]. The forthcoming experiments at FAIR, J- 
PARC and the low energy programs at RHIC would be infor- 
mative for these intermediate and high density states of matter, 
which may be realized in the neutron-star core. 

The QCD phase transition has two aspects; the chiral 
and deconfinement transitions. These transitions are asso- 
ciated with the spontaneous breaking of the global symme- 
try SU(N f ) L x SU(Nf) R in the chiral limit and the global 
Zjv e symmetry in the heavy quark mass limit, respectively, 
where the order parameters are the chiral condensate and the 
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Polyakov loop. Theoretical frameworks describing the dy- 
namics of these order parameters would be mandatory to un- 
derstand the QCD phase transition and phase diagram. The 
lattice QCD MC simulation is the most reliable and rigor- 
ous framework to investigate the non-perturbative aspects of 
QCD [5]. Recent MC simulations demonstrate that the two 
transitions occur simultaneously within the error bars at small 
or zero chemical potential @]. There is no clear explana- 
tion for the mechanism of coincidence of transitions, and 
some models suggest that the chiral and deconfinement tran- 
sitions can start to separate as quark chemical potential be- 
comes large. Unfortunately, MC simulations have the notori- 
ous sign problem at finite fj, [Q], and cannot answer whether 
two transitions can be separated or not. As a result, we need 
to invoke some approximations in QCD or effective models at 
finite chemi cal p otentials, such as the strong coupling lattice 
QCD S[il[il[ll[i3M|2iM[ll[MHl, Nambu-Jona- 
Lasinio model (NJL) 1321 and Polyakov loop extended NJL 
model (PNJL) l33T[34ll . and chiral random matrix model @]. 

Strong coupling lattice QCD (SC-LQCD) can provide a 
simple and lattice-based description for the chiral and decon- 
finement transitions at finite T and fi. This method is based on 
the strong coupling (1/g 2 ) expansion, where the link variables 
are integrated out analytically and it is possible to elucidate 
important aspects of the QCD phase transitions beyond a con- 
tamination of the sign problem. In the pure Yang-Mills the- 
ory, the color confinement in vacuum was shown in the strong 
coupling limit via the area law for the Wilson loop fioll . SC- 
LQCD of the pure Yang-Mills theory with higher order cor- 
rections ifTlll is recently extended to the deconfinement tran- 
sition study at finite T lfl2ll . At finite T, the deconfinement 
transition was qualitatively explained on the base of the mean- 
field treatment of the Polyakov loop in the leading order of 
the strong coupling expansion 11311 . and higher order correc- 
tions on the Polyakov loop action has been investigated re- 
cently lfl2tl . The deconfinement transition is also investigated 
by taking account of the Haar measure contributions to the ef- 
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fective potential 11411 . The behavior of the Polyakov loop is 
investigated also in the weak-coupling regime 1151 . 

We find some pioneering works on the chiral symmetry 
breaking in vacuum |36| and its restoration at finite T and/or 
p 11161 — TlSfl where one can find a sophisticated expression of 
the effective potential describing the chiral transition in the 
strong coupling limit Il6l - [l8ll . It is remarkable that a promis- 
ing phase diagram structure has been obtained even in the 
strong cou pling limit (SCL) analytically J 1^ 22| and in MC 
simulation 12311 . In our previous works I24ll25ll . we have de- 
veloped a formulation to evaluate the next-to-leading order 
(NLO, 0{l/g 2 )) Gil III Elk and the next-to-next-to-leading 
order (NNLO, (0(1/ a 4 )) i lHI^ of the strong coupling ex- 
pansion. Finite coupling effects appear as modifications of 
the quark mass, quark chemical potential, and the wave func- 
tion renormalization factor. The phase diagram evolution with 
j3 — 2N c /g 2 can be interpreted in terms of those modifica- 
tions. The SC-LQCD technique developed in these works is 
also applied to the graphene system in the condensed-matter 
physics l28ll . where the strong coupling gauge theory is appli- 
cable at low energies. 

In the works mentioned above, either of the two order pa- 
rameters is included, and the interplay between these order 
parameters is not discussed. One of the interesting develop- 
ments can be found in the works by Gocksch and Ogilvie 12911 
and by Ilgenfritz and Kripfganz 13011 . They developed a model 
including both of the chiral condensate (cr) and Polyakov loop 
(£). In this model, abbreviated as the GOIK model, the decon- 
finement is governed by the effective action for I derived in 
the leading order of the strong coupling expansion in the Pure 
Yang-Mills theory, and the the chiral transition is described by 
the strong coupling limit effective action for quarks l29l - [3lll . 
The coupling between a and I is found to appear naturally 
from the quark determinant. This coupling leads to the cor- 
relation of the chiral condensate and Polyakov loop, and the 
two transitions take place at similar temperatures 13111 . This 
understanding in the GOIK model led to the development 
of the Polyakov loop extended Nambu-Jona-Lasinio (PNJL) 
model il. 

From a view point of the rigorous strong coupling expan- 
sion, there is a problem to define the expansion order when we 
include both of the chiral condensate and Polyakov loop; the 
quark effective action has the strong coupling limit O(l/g ), 
while the leading order Polyakov loop action is proportional 
to l/g 2Nr , where N T is the temporal lattice size. We there- 
fore should regard the order of the strong coupling expan- 
sion separately in the quark and Polyakov loop effective ac- 
tion. Namely we start from the leading order in each of the 
quark and Polyakov loop action, O(l/g ) and (D(\/ g 2N ^), 
respectively, and extend the framework by including higher 
order terms of the strong coupling expansion in two direc- 
tions. As shown in Table U the effective action adopted in the 
GOIK model is the starting point, where both of the quark and 
Polyakov loop effective action terms are in the leading order. 
In this effective action, the QCD coupling constant appears 
only in the coefficient of the Polyakov loop action, and that 
coefficient J = (l/g 2 N c ) T is replaced with an expression 
J ~ exp(— na 2 /T) using the string tension k 129T - B1I1 . This 



is not fully consistent with the quark effective action where 
the coupling constant is assumed to be infinite. Therefore it 
would be more favorable to combine the quark effective ac- 
tion in SC-LQCD with finite coupling effects [nlElSl with 
the Polyakov loop action. Since the NLO and NNLO effects 
in SC-LQCD are represented as modifications in the form 
of the SCL effective action ll24i I25K we can combine them 
with the Polyakov loop effects in a similar way to that in the 
GOIK model. This Polyakov loop extended strong coupling 
lattice QCD (P-SC-LQCD) would be a promising framework 
to investigate the relation between the chiral and deconfine- 
ment transitions, and has some advantages over other effec- 
tive models. For example, P-SC-LQCD is directly based on 
the lattice QCD, and we can compare its results with those in 
MC simulations at zero chemical potential. 



In this paper, we develop a P-SC-LQCD framework by 
combining the leading order Polyakov loop effective action 
and NNLO quark effective action. We derive an analytic ex- 
pression of the effective potential in P-SC-LQCD at finite T 
and fx, and investigate the chiral and deconfinement phase 
transitions at p = 0. The whole phase diagram structure 
including finite p is soon reported elsewhere 1351 . Here we 
adopt one species of unrooted staggered fermion which cor- 
responds to four flavors in the continuum limit. The unrooted 
staggered fermion is suitable to developing an analytic for- 
mulation, since its simple structure enables us to obtain the 
effective potential analytically. The effective action is com- 
posed of the fermionic and pure gluonic parts. We adopt the 
fermionic effective action with finite coupling effects ifTvl |24| — 
l27ll . The leading order Polyakov loop action is obtained from 
N T plaquette configuration as shown in Fig. Q] We evaluate 
the Polyakov loop effects and derive the effective potential 
in two kinds of methods. One is the Haar measure method 
where we replace the Polyakov loop with its mean-field value. 
Instead of integrating out the temporal link variables, the log- 
arithm of the Haar measure is included in the effective poten- 
tial. The other is the Weiss mean-field method, which includes 
the fluctuation effects of the Polyakov loop. In this method, 
we bosonize the Polyakov loop action and carry out the tem- 
poral link integral. On the basis of these effective potentials, 
we investigate the chiral and deconfinement phase transitions 
simultaneously. We also compare the results of the two meth- 
ods. 



This paper is organized as follows. In Sec. [TTJ we derive 
the effective potential with Polyakov loop and finite coupling 
effects in both Haar measure and Weiss mean-field methods. 
In Sec.|nll we investigate the chiral and deconfinement phase 
transitions, and compare the results in the two methods. In 
Sec. lIVl we summarize our work and give a future perspective. 



Throughout this paper, we utilize the lattice unit a = 1, 
and physical values are represented as dimensionless values 
normalized by the lattice spacing a. 
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TABLE I: Classification of SC-LQCD at finite T and fj,, The unrooted staggered fermion is utilized in the quark effective action except 
for fl3l . where the Wilson fermion is included in the hopping parameter expansion. In Ref. lfl3h . they treat the chiral spin model, which is 
essentially the same as that of the Polyakov loop action. 





Polyakov loop action 


Chiral quark action 


w/o Polyakov loop 


Leading Order 


Higher orders 


w/o quarks 




Kogut, Snow, Stone 1 13] 
Polonyi, K. Szlachanyi; Gross; 
Bartholomew et al. fl4ll 


Langelage, Munster, Philipsen f 121 


LO (SCL, O(l/g )) 


"SCL-LQCD" 
Refs. (USB 


"GOIK model" 
Gocksch, Ogilvie I29fl 
Ilgenfritz, Kripfganz 13011 
Fukushima ll3lll 


— 


NLO (0(l/g 2 )) 


"NLO SC-LQCD" 
Faldt, Petersson fl7ll 
Bilic et al. |26[1 
Miura, Nakano, Ohnishi, 
Kawamoto l24ll 


"P-SC-LQCD" 
Nakano, Miura, Ohnishi 
(present work) 
Miura, Nakano, Ohnishi (3^1 




NNLO (0(l/g 4 )) 


"NNLO SC-LQCD" 
Jolicoeur et al. j27h 
Nakano, Miura, Ohnishi |2^I 





H. EFFECTIVE POTENTIAL 

A. Lattice QCD action 

In the lattice QCD, the partition function and action with 
one species of unrooted staggered fermion for color SU (N c ) 
in the Euclidean spacetime are given as, 

Zlqcd = J V[x, X, U v ] e~ s ^ CD , (1) 
Slqcd =Sf + Sg , (2) 
1 d 

S f =r [v P ,xXxU p , x Xx+p ~ V P ,lXx+pUl x Xx] 



X p=Q 



u P 



(3) 
(4) 



where x(x)> TO o> U px , and Up denote the quark (antiquark) 
field, the bare quark mass, link variable, and plaquette, respec- 
tively. In the staggered fermion, the spinor index does not ap- 
pear explicitly and its structure is expressed in the staggered 
phase factor n P>x = (vo,x,r]j, x ) = (e", (-If o+-+xj-i) 
13811 . The quark chemical potential \i on the lattice is intro- 
duced as a weight of the temporal hopping in the staggered 
factor ll39ll . The above lattice QCD action has a global sym- 
metry U(1)l X U(1)r in the chiral limit (mo — > 0), and a stag- 
gered chiral transformation is given as y z — > e x \x where 
e x = (-l)*o+~+x d i s re i a ted to 75 iHlgl. 

Throughout this paper, we consider the case of color 
SU(N C = 3) in 3+1 dimension (d = 3) spacetime. Tempo- 
ral and spatial lattice sizes are denoted as N T and L, respec- 



tively. While T = 1/N T takes discrete values, we consider 
T as a continuous valued temperature. We take account of 
finite T effects by imposing periodic and antiperiodic bound- 
ary conditions on link variables and quark fields, respectively. 
We take the static and diagonalized gauge (called Polyakov 
gauge) for temporal link variables with respect for the peri- 
odicity iflql . In these setups, we derive the effective poten- 
tial J-" off = — log[ZLQco] / {N T L d ) on the basis of the strong 
coupling expansion with Polyakov loop effects. 



B. Effective action with Polyakov loop effects 

To treat the chiral and deconfinement phase transitions si- 
multaneously, we evaluate the leading order Polyakov loop 
effects in the pure Yang-Mills sector and NNLO effects in the 
fermionic sector. In a finite T treatment of SC-LQCD, we 
first derive an effective action from the spatial link (Uj) inte- 
grals, and the temporal link (Uq) integral is evaluated later to 
consider the thermal effects of quarks |[T7i[l8ll . 

We consider the leading order Polyakov loop effective ac- 
tion Dg ff of the strong coupling expansion in the pure Yang- 
Mills sector, which is obtained from the spatial link integral 
of the N T plaquettes sequentially connected in the temporal 
direction as shown in Fig.Q] 



'off 



1 



N- 



xj>0 



P.P. 



h.c 



■)■ 



(5) 



where P x = tr c Y[ T Uq(t, x)/N c represents the Polyakov 
loop. The factor l/g 2 N c arises from the spatial link integral. 
This effective action shows the nearest-neighbor interaction 
of the Polyakov loop. While we cannot reach the Stefan- 
Boltzmann limit in the high temperature limit with lfl2ll . 
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it is found to give results qualitatively consistent with the 
MC simulations on the first order phase transition properties 
such as the deconfinement temperature and the discontinu- 
ity of the Polyakov loop in the mean-field treatment of SC- 
LQCD Il3l Il9ll . then essential dynamics around the decon- 
finement transition regime may be described. 



x x + j 



D:(Mn) 

spatial ~ 

P P 

link -fx r x+j 

integral 



FIG. 1: Leading order of the Polyakov loop effects in the strong 
coupling expansion. The squares in the left and loops in the right 
represent the temporal plaquettes and the Polyakov loops, respec- 
tively. 



As for the effective action for quarks, we utilize the NNLO 
effective action ll25ll derived from one- and two-plaquette 
configurations, which correspond to 0(l/g 2 ) and 0(l/g 4 ), 
respectively, in addition to the strong coupling limit terms 
(0(1/ g )). We consider only the leading order terms of the 
1/d expansion, 0(l/<i o ), which corresponds to taking mini- 
mum quark number configurations for a given plaquette ge- 
ometry ll-ioTl . We briefly explain the NNLO effective action in 
AppendixlAl 

The effective action including both of the quark and 
Polyakov loop contributions is represented as 



5cir 



+ S, 



(P) 



S'gg represents the NNLO fermionic effective action, 



G^(m q ; (M,T) = - [e»U 0iX 5 x+ ^ y e""^ 8 X _^ 



(6) 



(7) 



(8) 



Finite coupling effects result in the modifications of the wave 
function renormalization factor Z x , quark mass rh q , and ef- 

( x) 

fective chemical potential jl. S cS represents the auxiliary 
field part of the effective action described in the AppendixlAl 

From a naive counting of the strong coupling expansion 
order, the leading order Polyakov loop action composed of 
the plaquettes is in the higher order (0(l/g 2NT )), compared 
with the NNLO terms {0(1/ g 4 )) stem from the quark sec- 
tor. As discussed in Introduction, we set the starting point of 
chiral and deconfinement dynamics as the effective action of 
0(1/ g°) and 0(l/g 2Nr ) in the quark and Polyakov loop sec- 
tors, respectively. The effective action in Eq. © corresponds 
to an extension in the quark sector from this starting point. 

From these action terms, an approximate QCD partition 
function is obtained, and the effective potential is defined by 



the logarithm of the partition function as 
1 



cir 



los 



N T L d 

=^;H,T) + U g {i,l)+J^ ) ^) 



(9) 



Here J" q ($; fj,, T), U g (£, I), and ) ($) represent the quark 
free energy, pure gluonic potential, and the effective potential 
including only the auxiliary fields, respectively. We obtain 
Jq(<l>; fx, T) and U g (£, I) by evaluating the Grassmann (x, x) 
integral and the temporal link (Uq) integral, 

J VU e- s ^Bet [Z x G _1 (m,,;/S,T)] 

= Y[ e N e (lo sZx+ E q )/T f d^e-Sg* 
x J 

xdet c [(l+W e-^- W/T ) (l+W t e-( B «+«/ T )] , 

(10) 

where Uq (x) = n r ^o(x,T). Note that the determinant is 
for the spacetime and color in the first line of Eq. ( [Tol l, and 
for the color in the second line of Eq. ( flOj) . In Eq. ( [Tol l, 
we carry out these integral as follows ll6lll7il2Tll . First, we 
perform the Fourier transformation in the temporal direction, 
and obtain the product in the frequency using the Grassmann 
integral over the quark fields. Second, we evaluate the product 
in the frequency by the Matsubara method. Finally, in this 
paper, we evaluate the temporal link integral in two kinds of 
schemes, the Haar measure and Weiss mean-field methods. 
We explain these methods in the next two subsections. 



C. Haar measure method 

We shall now derive the effective potential with Polyakov 
loop effects in the Haar measure method (H-method). In 
the H-method, we replace the Polyakov loop with a mean- 
field value and the Haar measure is taken into account in the 
Polyakov loop potential instead of carrying out the temporal 
link integral. The contribution to the effective action is, 



S. 



(P) 

efi 



2f3 p L d ££ . 



(11) 



where (3 p = (l/g 2 N c ) N - N 2 d_, I = (P x ) and I = (P x ). We 
assume the mean fields £ and £ are constant and isotropic. 

The temporal link integral is represented as the Haar mea- 
sure in the Polyakov gauge, 



dU = / d£d£ ■ 27 



1 - 6££ + 4 (£ 3 + F) -3(£i)'' 

(12) 

The Polyakov gauge is a static and diagonalized gauge for 
temporal link variables lfl6ll . then Uq is given as follows, 

W (x) =Y[u (x,t) =dmg c (e te \e^,e^) . (13) 
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This Haar measure shows the Jacobian in the transformation 
from the temporal link variables (Uq) to the Polyakov loop 
(£, £). Therefore, J" q ($; (j,, T) and U g {£, £) are given as, 

Jq = -N c E q - T log R(E q - fi, N c £, Njt) 

- T log R(E q + fl, Nj, N c £) - N c log Z x , (14) 
R(x, L, L) = 1 + Le~ x l T + Le~ 2x / T + e~ 3x / T , (15) 

U g = -2T/3 p l£ - Tlog [l - 6£l + 4 (£ 3 + F) - 3 (££)' 

(16) 

where f3 p = (1/ g 2 N c ) 1 / T N 2 d. Here we have replaced the 
N T with 1/T. Note that we have omitted irrelevant constants. 
In Eq. (1141) . the first line represents the vacuum and quark free 
energy and the second line represents the antiquark free en- 
ergy and the contribution of the wave function renormalization 
factor. The quark free energy includes one- and two-quark 
excitations (e~( E ~^/ T , e~ 2i > E ~' :L ^ T ). In the confined phase 
(£ ~ I ~ 0), the one- and two- quark excitations are sup- 
pressed and only the color-singlet state contributions remain. 
The pure gluonic potential U g does not include the fluctuation 
of the Polyakov loop since we treat the Polyakov loop as the 
mean field without the temporal link integral. In the PNJL 
model, this pure gluonic potential is incorporated to express 
the properties of the deconfinement phase transition l33ll . 



one-quark excitation two-quark excitation 



(EH 




one-antiquark excitation 




temporal baryon hopping 



I 



FIG. 2: Quark excitation on the lattice. We show the one- and two- 
quark excitation in the first line, while the temporal baryon hopping 
in the third line. 



transformation is a procedure to bosonize the product of dif- 
ferent types of composites 12411 . The Polyakov loop action in 
Eq. ^ is linearized as, 



S. 



(P) 

eff 



1 



9 2 N C 



N 2 2(M-P X £-£P X 

x,j>0 

2p p L d l£ - 2/3 p Y, ( P * £ + lP *) > 



(18) 



where f3 p ,P x , and P x are defined in Subsec. HI B I £ and 
£ represent the auxiliary fields for the Polyakov loop, {£ = 
(P x },£ = (Px))- In the last line of Eq. HSi . we assume 
constant and isotropic values for auxiliary fields £ and I. 

After the Grassmann integral, we carry out the tem- 
poral link integral in Eq. ( fTOb explicitly in the Polyakov 
gauge. Since includes the Polyakov loop (P X ,P X ) = 

we have to include S^P in the tempo- 
ral link integral in the W-method. The relevant part of the 
temporal link integral in Eq. ( fTOb is given as, 



Z P = dU det c 



2 cosh{E q /T) + U^ IT + z4 e ~ A/T 



x exp 



p ■ 



(19) 
(20) 



We have defined 77 = 2/3 p £/N c and 77 = 2/3 p £/N c . In Ap- 
pendix [B] we explain the detail of the temporal link integral. 
Eventually, we derive the effective potential, 



T q = - N c E q - T log(Z' P /L ) - N c log Z x 



(21) 



Z' P /Lo =R(E q - ft, L lt Li) R(E q + £, LxM) (22a) 



,-2EjT 



1 + e- 2E "/ T 



{2 + 12-1,1,) (l + e- 2E ^ T 
(2L X + L 3 - L 2 ) e-^/ T 

(2Li + U - L 2 ) e^ E+ ^' T 



(22b) 

(22c) 
(22d) 

(22e) 



D. Weiss mean-field method 

Now we shall evaluate the Polyakov loop effects in the 
Weiss mean-field method (W-method). The W-method in- 
cludes some part of the fluctuation effects of the Polyakov 
loop. We first bosonize the Polyakov loop action by using the 
Extended Hubbard-Stratonovich (EHS) transformation 12411 . 

e aAB w cxp [-a(4>i> -Aip- $B)] , (17) 

where a denotes a positive constant, and two auxiliary fields 
are introduced simultaneously, {^>,ip) = ((A), (B))). EHS 



U 9 {£, I) =2Tp p M - T log L Q . (23) 

Lq, Li, Li, L2, L3, L3 are the functions of 77 and fj. These 
functions include the modified Bessel function /„ (2^/7777) and 
the ratio of £ and I (<j> = §log(7?/?7) = ±log(£/Y)). The 
explicit expressions of them are shown in the AppendixlBl 

Compared with the effective potential in the H-method 
where the quark contributions are represented as the vacuum, 
quark and antiquark parts , Eq. ( fT~4T >. the effective poten- 
tial is more complicated. While Eq. d22ab represents quark 
and antiquark parts (including the one- and two-quark exci- 
tations, and temporal baryon hopping) as H-method, we re- 
place £{£) with Px(-Px) in Fig. [2] . The combination of 
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the Bessel function accompany with the Boltzmann factors 
in W-method. Other terms represent the temporal meson hop- 
ping (Eq. (122bl i). one- or two- quark and one- or two- anti- 
quark excitations (Eq. (122cl i), and one-antiquark (one-quark) 
and two-quark (two-antiquark) excitations (Eqs. (122dl i,( 122el )). 
as shown in Fig. [3] We can understand these physical in- 
terpretations by considering the Boltzmann factors. We ex- 
plain the one-quark and one-antiquark excitations {e~ 2Eq / T ) 
in Eq. ( I22cb as an example. We can decompose Eq. (122al i. 
and find that the contribution of (e~ 2E i/ T ) comes from the 
multiplication of one-quark (e~( Eq ~^/ T ) and one-antiquark 
^ e -(E q +p,)/T^ exc itations. These contributions of Eqs. ( 122bl >- 
(I22eb in W-method can appear because of the temporal link in- 
tegral and not in H-method or the PNJL model. In the confined 
phase (£, I ~ 0), the quark excitations are suppressed since 
Li ^ Li ~ L3 ~ L-3 ^ £2 + 2^0, and the color-singlet 
state is dominant. Therefore, Eqs. ( I22ab and (I22bb remain fi- 
nite in the confined phase, while the other terms Eq. ( 122c| >- 
( |22et vanish. 

Note that the pure gluonic potential U g includes the depen- 
dence on £/£ explicitly (i.e. the dependence on /i). When the 
quark chemical potential is zero (/i = 0), the Polyakov loop 
for quarks and antiquarks are the same (£ = £). In compari- 
son, when the quark chemical potential is finite (/1 7^ 0), the 
Polyakov loop for antiquarks is generally different from that 
for quarks (£~f £) & 

In Ref. IU3I1 . the pure gluonic potential is derived at fj, = 
without quarks. As shown later, we confirmed that our results 
are consistent with their results, for example, in the tempera- 
ture of the deconfinement phase transition without quarks. 

(b) (d) 
temporal meson hopping two-quark one-antiquark 
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(c) (c) 
one-quark one-antiquark two-quark two-antiquark 
0/^\« »'^\Q 00« 

• »_/o o >. 

FIG. 3: (b)temporal meson hopping (c)one-quark and one- 
antiquark excitations, two-quark and two-antiquark excitations 
(d)one-antiquark and two-quark excitations 



III. CHIRAL AND DECONFINEMENT PHASE 
TRANSITION 

We shall now discuss the phase transitions in P-SC-LQCD. 
We impose the stationary condition for the effective potential 
with respect to the auxiliary fields, and obtain thermodynam- 
ical quantities. In this article, we investigate the chiral and 
deconfinement transitions at N c = 3 and /1 = 0. First, we 



discuss the deconfinement transition in pure Yang-Mills the- 
ory, and next we show the results in the W-method. Finally 
we compare the results in the W- and H-methods. 



A. Deconfinement transition in pure Yang-Mills theory 



P=2N c /g 2 =4.0, T/T d c =0.9,1.0,l.l 



0.5 



-0.5 



Weiss 
Haar 

T dc (Weiss) =0.7491 




Haar 

T dc =0.6127 




0.5 

£ 



FIG. 4: Polyakov loop potential in the W- (upper panel) and H- 
(lower panel) methods. The filled circles show the equilibrium 
points. We show the results in the lattice unit. 

Before discussing the results with quarks, we first examine 
the deconfinement transition in the pure Yang-Mills theory. In 
Fig. |U we compare the Polyakov loop potentials U g in the W- 
and H-methods, Eqs. d23l ) and ( fToT i. respectively. We show the 
results at T = (0.9 - l.l)T djC at p = 4 without quark effects. 
Both of the potentials have qualitatively the same behavior 
and show the first-order transitions at T = Td. c . We find some 
quantitative differences in these two potentials. First, we note 
that in H-method there is a wall at £ = 1 where U g diverges. 
The range of £ is limited in H-method to be less than one as 
can be found from Eq. (116t . In the W-method, the integral 
over Uq smoothens the H-method potential, 



(U)/t = / de i e -/3p(e-e.')(£-£') e -u% '(e ',e')/T 



(24) 



and the effective potential becomes a milder function of £. As 
shown in the dotted line in the upper panel of Fig. |4j is 
smaller (larger) for large (small) £. Since the transition tem- 
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TABLE II: Correspondence of the parameters between our results 
and Ref. (7J]. Note that we replace /3 in Ref. {H with /3 (K) to 
distinguish from j3 — 2N c /g 2 . 



Our results Kogut et al. [TJ] (d = 3) 

£ a/2d 
2f3 p /N 2 c 2d/3 ( - K) /(2N c ) 



perature is more sensitive to the potential at small £, the T^ c 
in the W-method is higher than that in the H-method. 

Our treatment in W-method in the pure Yang-Mills theory 
is the same as the chiral spin model investigated in Ref. Il3ll . 
where the following effective action is adopted, 



'oft 



/L d = 



Ad 



a(P + P) 



(25) 



When the parameter /3( K ' and the variable a are rewritten ac- 
cording to Table [Tj] this action is found to be equivalent to 
S^P in Eq. ( fl"8l ). In Ref. iQ, the first-order phase transition 
is found to take place at the critical coupling /?, 



(K) 



0.81 



in the mean-field treatment. By using the relation (3^ = 
2(3 p /N c d, this critical coupling corresponds to T^e ~ 0.75 at 
(3 = 4. This value is consistent with our result Td c ~ 0.7491 
in the upper panel of Fig. [4] 



B. Chiral condensate and Polyakov loop 



P = 4.0 




[3 = 4.0 




FIG. 5: Temperature and quark mass dependence of the chiral con- 
densate (upper panel) and the Polyakov loop (lower panel) in the 
W-method. We show both of the results in the lattice unit. 



We shall now discuss the chiral and deconfinement transi- 
tions in the W-method at finite T and /j, = 0, via the T depen- 
dence of the chiral condensate a = — (xx) an d me Polyakov 
loop (£, £) = ((-P x ), (Px))- In Fig. [3] we show the quark-mass 
(too) and T dependence of a and I at = 0. We adopt the 
coupling /3 = 4 as a typical example. As expected, a vanishes 
sharply at T c at too = showing the second order transition 
in the chiral limit, and £ starts to show the first order transition 
at mo ~ 4.0. 

In Fig. [6] we show a and I as functions of T at /? = 4 and 
mo = 0.05. Since the quark mass is small but finite, both of 
a and £ show crossover transition towards the chiral restored 
(cr — > 0) and deconfined (£ — > 1) matter, i.e. a quark-gluon 
plasma. 

From Fig.|6j the chiral and deconfinement transitions seem 
to take place almost at the same temperatures. In the upper 
panel of Fig. [7] we show —da/dT and d£/dT as functions of 
T. Both of -da/dT and d£/dT are peaked at around T ~ 
0.55. We here define the chiral T x c and deconfinement Td , c 
transition temperatures as the peak temperatures of —da/dT 
and d£/dT. In the lower panel of Fig. [7] we display T XjC and 
Td.c as functions of (3. 

Polyakov loop is found to suppress the chiral conden- 
sate and to reduce the chiral transition temperature in P-SC- 
LQCD. In Figs. [6] and [8] we compare the chiral condensate 
a and the chiral transition temperature T x , c with and without 



TABLE III: Difference between the SC-LQCD with/without the 
Polyakov loop (PL) effects, the PNJL, and NJL model. 



low T high T 

SC-LQCD w/o PL confined 
NJL deconfined 
P-SC-LQCD and PNJL confined deconfined 



Polyakov loop effects. We find that both a and T XjC in P-SC- 
LQCD are smaller than those without Polyakov loop effects; 
the chiral condensate becomes smaller when the Polyakov 
loop takes a finite value, and the chiral symmetry restoration 
takes place at lower T by the Polyakov loop effects. In SC- 
LQCD without Polyakov loop effects, we find the contribution 
only from color-singlet states as discussed before, then we im- 
plicitly assume that the quarks are confined at any T. In the 
W-method, we have one- and two-quark contribution as well 
when the Polyakov loop takes a finite value. Quark excita- 
tion generally breaks the chiral condensates, then it promotes 
the chiral symmetry to be restored at lower T. The similar 
behavior is found also in the H-method. 

This behavior is different from that in the PNJL model 
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P=2N c /g 2 =4.0, n=Q, m =0.05 
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P=2N c /g 2 =4.0, n=0, m =0.05 
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FIG. 6: Chiral condensate and Polyakov loop in P-SC-LQCD (solid 
lines), and chiral condensate in SC-LQCD without the Polyakov loop 
effects (dashed line) as functions of T at /j, — in the W-method. We 
show the results in the lattice unit. 



033U34H . In PNJL, the Polyakov loop pushes up the chiral tran- 
sition temperature from that in the NJL model. This difference 
can be understood as follows. Quarks are deconfined in NJL, 
then the role of the Polyakov loop is to confine quarks at low 
T in PNJL. By comparison, only color-singlet hadronic states 
contribute to the effective potential and quarks are confined in 
SC-LQCD w/o PL, then the role of the Polyakov loop is to 
deconfine quarks at high T in P-SC-LQCD. As a result, con- 
fined and deconfined phases appear at low and high T region 
in both of PNJL and P-SC-LQCD, and the upward and down- 
ward shifts of T c in PNJL and P-SC-LQCD are understood as 
the directions to enlarge the region which is promoted by the 
Polyakov loop. As a result, both of PNJL and P-SC-LQCD 
give a consistent picture of the chiral and deconfinement tran- 
sitions, although the confinement properties of quark matter 
are different in NJL and SC-LQCD without Polyakov loop ef- 
fects as shown in Table ITTT1 

Note that the confinement picture of quark matter in PNJL 
model and P-SC-LQCD is not completely the same. In P- 
SC-LQCD, colors of quarks in a hadron are confined locally 
when Polyakov loop are zero. By comparison the inter-quark 
distance in a hadron is not necessarily small in PNJL, since 
the momentum of quarks are specified. 

In our previous work on NLO H and NNLO J3 SC- 
LQCD, the critical temperature is calculated to be larger than 
MC results, and NNLO effects on T X:C are found to be small. 
In Figs. [8] and [9] we show T XtC in several treatments of P- 
SC-LQCD in comparison with the MC results. As already 
discussed, the chiral transition temperature T x . c is reduced by 
the Polyakov loop effects and roughly explains the MC re- 
sults, especially in the region j3 < 4. By comparison, the 
NNLO effects on T X}C at /x = are small also in the case with 
Polyakov loop effects. This observation implies that introduc- 
ing the Polyakov loop, the deconfinement order parameter, is 
essential to explain the QCD phase transition temperature. In 
order to explain further reduction of T x , c at f3 > 4, we may 
need to include higher order terms of 1/g 2 in the Polyakov 
loop effective action. 
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FIG. 7: Upper panel: Temperature dependence of —do/dT (dashed 
line) and dl/dT (solid line). Lower panel: /3 dependence of the 
critical temperatures for the chiral (dashed line) and deconfinement 
(solid line) phase transitions. All of them are obtained in the W- 
method. We show both of the results in the lattice unit. 



Comparison of Weiss mean-field and Haar measure 
methods 



The Haar measure method is useful to include the confine- 
ment effects in chiral effective models, and it has been widely 
applied for example in the PNJL flU and PQM JH models. 
In this subsection, we discuss the results of the Haar measure 
method (H-method), and compare the results of two methods. 

First we compare the chiral condensate and the Polyakov 
loop in the W- and H-methods. In Fig. [TO] we show the T de- 
pendence of a and I in the two methods. While qualitative be- 
haviors are the same in both of the methods, the chiral conden- 
sate (Polyakov loop) in the W-method is larger (smaller) than 
those in the H-method, suggesting that the W-method exhibits 
a little larger T c than the H-method. In Fig. [8] we compare 
the chiral transition temperature T x c in the two methods. The 
transition temperature is defined from the peaks of —da jdT 
and d£/dT as shown in Fig.fTTI 

Physically, the temporal link integral in the Weiss mean- 
field method favors the color-singlet states and therefore sup- 
press the quark excitation. Also, as another aspect, the smaller 
Polyakov-quark coupling in the W-method is responsible in 
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FIG. 8: Comparison of P-SC-LQCD (solid line) and SC-LQCD without the Polyakov loop (dashed line) in the lattice unit. We show the results 
of Weiss mean-filed method (bold solid line) and Haar measure method (thin solid line) in P-SC-LQCD at mo = 0.0. The triangles represent 
the results of the critical temperature (T c , M= o, open triangle) and the critical coupling (/3 C , filled triangles) obtained in Monte-Carlo simulations 
with one species of unrooted staggered fermion: From the left, T Ct)1= o in the SCL with monomer-dimer-polymer (MDP) simulations l23l . /3 C 
at (N T ,m ) = (2, 0.025) 0, (2, 0.05) 0, (4, 0.0) 0], (4, 0.05) H0, and (8, 0.0) 0]. 
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FIG. 9: Comparison of the chiral transition temperature T x , c in P- 
SC-LQCD with NNLO (solid line) and NLO (dashed line) finite cou- 
pling corrections for quarks. The triangles represent the same MC 
results in Fig. [8] We show the results in the lattice unit. 



part for the larger T c . In the W-method, we carry out the tem- 
poral link integral of the effective action in which fermions 
couple with the temporal link variables. This "averaging" 
procedure as shown Eq. d24l i smears the Polyakov-quark cou- 
pling. In the upper panel of Fig.Q~2] we show the comparison 
of the Polyakov-quark couplings, Li(= L{) and N c l in the 
W- and H-methods, respectively, defined as the factor accom- 
panied with the Boltzmann factor. At \i = 0, we find that 
the stationary value of coupling is smaller in the W-method, 
Li(= Li) < N c £, at any temperature as shown in Fig. [12] 
These observations confirm the smaller Polyakov-quark cou- 
pling. Similar behavior is also found in the NJL model with 
the Polyakov integration 14711 . 



IV. CONCLUDING REMARKS 

In this paper, we have derived an analytic expression of 
the effective potential at finite temperature and chemical po- 
tential in the strong coupling lattice QCD with the Polyakov 
loop effects using one species of unrooted staggered fermion. 
We have discussed the chiral and deconfinement transitions at 
[i = 0, especially the Polyakov loop effects on these transi- 
tions. We have considered NNLO effects in the strong cou- 
pling expansion and the leading order of 1/d expansion in 
the quark effective action. We have adopted the leading order 
Polyakov loop action, 0(l/g 2NT ), and the temporal link (Uq) 
integral is evaluated in two methods. One is the Haar measure 
method (H-method), where we replace the Polyakov loop with 
its constant mean-field without the Uq integral. The decon- 
finement dynamics is taken account of via the Haar measure. 
Another is the Weiss mean-field method (W-method), where 
we bosonize the effective action using the Extended Hubbard- 
Stratonovich transformation and carry out the temporal link 
integral explicitly. The fluctuation effects appear as the com- 
bination of the modified Bessel function and the difference of 
£ and I in the effective potential. We have found that one- and 
two-quark excitations are allowed in both methods, when the 
Polyakov loop takes a finite value at high T. 

We find that the Polyakov loop reduces the chiral transition 
temperature, and the obtained transition temperatures roughly 
explain the MC results in the region (3 < 4. The chiral and de- 
confinement transitions are found to take place at similar tem- 
peratures, but the two transition temperatures are slightly dif- 
ferent when we define them from the susceptibilities. We have 
compared the results in the W- and H-methods. These meth- 
ods exhibit qualitatively the same results, but the Polyakov 
loop effects are found to be weaker in the W-method. The 



10 



0.5 




1 



0.5 



p=2N c /g 2 =4.0, n=Q, m =0.05 



Weiss 
Haar 



(3=2N c /g 2 =4.0 



0.5 
T 



FIG. 10: Chiral condensate and Polyakov loop as functions of T at 
/i = in the W- and H-methods. The notation is the same as Fig. [6] 
We show the results in the lattice unit. 



|3=2N c /g z =4.0, n=Q, m =0.05 



-d(o/N c )/dT 




FIG. 11: Temperature dependence of —da/dT (dashed line) and 
di/dT (solid line) in the Haar measure method (H-method). The 
notation is the same as Fig. [7] We show the results in the lattice unit. 



transition temperature is slightly higher and the Polyakov- 
quark coupling is smaller in the W-method. This is because 
the temporal link integral in the W-method favors the color- 
singlet states and leads to the suppression of the quark excita- 
tion. 

There are several points to be studied further. Firstly, we 
should investigate the QCD phase diagram in finite T and \i on 
the basis of the effective potentials derived in this work. Sec- 
ondly, it would be a promising direction to extend the present 




0.5 
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P=2N c /g 2 =4.0, 11=0, m =0.05 




FIG. 12: Upper panel: The Polyakov-quark coupling L\ and N c £ as 
functions of I at /i = (T=0.5,..,1.0) (solid line). The dots represent 
the value of L\ under the stationary condition (mo = 0.05). Lower 
panel: Stationary value of the Polyakov-quark coupling L\ and N c l 
as functions of the normalized temperature T/T c . We show both of 
the results in the lattice unit. 



P-SC-LQCD framework by including higher order terms in 
1/g 2 and 1 jd in the quark and Polyakov loop effective action. 
Thirdly, it is interesting to apply the present effective potential 
to the chiral effective model, for example to the PNJL model, 
in order to include the fluctuation of the Polyakov loop. We 
replace the Polyakov loop potential with that in Eq. d23l and 
the factor accompanied with the Boltzmann factor N c £ (N c £) 
with L\ (Li) in Eq. d2Tb . These studies will be helpful to 
shed light on the QCD phase transitions at finite T and yu. 
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Appendix A: NNLO effective action 

We here briefly introduce the NNLO SC-LQCD effective 
action l25ll . which we utilize as the quark part of the P-SC- 
LQCD. The NNLO effective action is derived from one- and 
two-plaquette configurations, which correspond to 0(1/ g 2 ) 
and 0(l/g 4 ), respectively, in addition to the strong coupling 
limit terms (0(1/ g )). We consider only the leading order 
terms of the 1/d expansion, O(l/d ). The spatial link inte- 
gral leaves a sum over spatial directions, and the energy per 
bond is assumed to be proportional to 1/d in order to keep the 
action finite at large spatial dimension in the 1 /d expansion 
lf4Qh . Since the quark field scales as d -1 / 4 , the leading order 
1/d terms correspond to the minimum quark number config- 
urations for a given pla quet te geometry. The NNLO effective 
action is represented as 12511 



o(NNLO) _ q (F) 
^eff —°efi 



s. 



iff j 



(Al) 



represents the fermionic effective action, shown in 
Eq. 0. 



D eff ~ 



\ E [Z.V+ - Z+V-] + J2 m' q M x 

X 

E \ [ e ~ s " v * + - e ^ r ] + E ™* M * 



Xv 



(A2) 



V+ =x x e»Uo, xXx+ s , V- = X x+6 e-^Ul xXx , (A3) 
M x =XxXx ■ (A4) 

Finite coupling effects result in the modifications of the wave 
function renormalization factor Z x = J Z+Z-, quark mass 
fh q = m' q /Z x , and effective chemical potential fl = fi — 5/i — 

/j, — log y/ Z + 1 Z- . Z± and m' q are defined as, 

m' q =b' a <r + m - l3 TT (i> TT + ip TT ) , (A5) 
Z+ =1 + /3'Jj T + Ap TT m'$ TT , (A6) 

Z_ =1 + (3' T ljj T + 4/3 TT TO^rr ■ (A7) 

Sgff "* represents the auxiliary field part of the effective action, 

=S$/{N T L d ) 

+Ptt4>tt^tt + fiss^sslpss + T^rsVVsVVs , (A8) 

Auxiliary fields (a, {ipK t4>k] K = t, s,tt,ts, ss}) in Eq. 
d A8t are introduced as mean fields for the mesonic compos- 
ites. We show the relation between auxiliary fields and the 
mesonic elements in Table [IV] The coefficients of the auxil- 



TABLE IV: The auxiliary fields in the NNLO effective poten- 
tial. Wx and WZ are represented as the mesonic composite 
which connects the quark fields in the next-to-nearest neighbor- 
ing temporal sites (W+ = x^Uo, x U 0x+i} x x+2 6 and W~ = 
Xx+2d e ~ 2ti U^ x+ qUq xX x)- We omit the indices x (x in the case 
of £, £) in the mesonic composites. 



Auxiliary Fields Mean Fields 



a 

(4>ss,4>ss) 
(VV,Vv) 

(1,1) 



—M 
(W+,W~) 
(MM, MMMM) 
(-V+V- ,2MM) 

(-v+,v-) 

(MM, MM) 



iary field in Eq. ( IA8b are defined as 

b' a = b„ + 2[j3 ss ^ ss +l3 TS i} TS +l3' s (>ip s +i> s )] , (A9) 

ft = Pr + Mrs , P' s = P S + 20„^„ , (A10) 



d d f 1 

a _ 2N C ' Pt ~ N 2 g 2 { + 2g 2 



(All) 



16Njg 4 



2N^ 



Appendix B: The temporal link integral in the Weiss mean-field 
method 



The part of the temporal link integral in Eq. ( fTOb is given 

as, 

Z P = J dU cxp [2P P (P X £ + IP X )] e N ^' T 
x dct c { [l +W e- ( ^-' i)/T ] [l +W V^ + «/ T ] } 

[dU ] Jj exp ( V e~ ie " + fje 1 ' 



a=l 



X + Ye 



i0„ 



(Bl) 



We define 7/ = 2fi p £/N c , fj = 2/3 p £/N c , X = 2cosh(N T E q ) 
and Y = exp(N T fi). 

This type of temporal link integral is written for N c = 3 
aslEl, 



dU ]]Jf{9 a )= E NJdetMn+i.j 

a—1 n— — 00 

N c \ E [Ml + M 2 n+1 M n _ 2 + M^_]M n+ 2 

n 

-M n (M n+2 M n _ 2 + 2M n+1 M„_i)] , (B2) 
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where 



M n = J ^/(0)exp[-m0 



X + Ye 



id 



Y 



(B3) 
(B4) 



where z = 2i^/ffr] and 4> = h ^ogiv/v)- By using to the re- 
lation e lzsin0 = Jii z )e iW > where M z ) is the Bessel 
function and I, z are an integer and a complex number, respec- 
tively, f(9) is found to be, 



f{6) = J2 



Jl(6-Tr/2-i4 



i=-c 



-10 



X + Ye ib 



Y 



(B5) 



We substitute Eq. ( |B5t into Eq. ( |B3t and carry out the integral 
over 0. The matrix element M n is calculated to be 



I n (y)X + I n ^(y)e^Y + I n+1 (y) 



Y 



(B6) 



Note that I v {y) = e"" 7 ™/ 2 J v (iy) is the modified Bessel 
function, where y = 2^/rfq. By substituting Eq. ( IB6b into 
Eq. ( IB2b . the temporal integral is found to be (up to a constant 

(NM 

Z P =D ± X 3 + D 2 X + D 3 Y 3 + D^ 3 Y- 3 

+X 2 (D 4 Y + D^Y- 1 ) + X(D 5 Y 2 + D^ 5 Y~ 2 ) 
+D 6 Y + D_ 6 Y-\ (B7) 

where Y = Fe -0 . Coefficients A (i = 1,2, ±3, • • • ,±6) 
are defined as, 



In the numerical calculation, we have confirmed that the 
sums in A converge. A.i (i = 1, 2, ±3, • ■ ■ , ±6) are rep- 
resented as combinations of the modified Bessel functions. 
From Eqs. (IB 1 U and dB13b . D±3 and D±$ are represented 

as D± 3 = e ±A ^A and D± 5 = e ±iV ^D T4 . We find 
A' (Vi v) = Di(fj, rf) from the relation £>-„.;' = A,i (£' = i 
for i = 1,2 and z' = — i for i = 3, • • • , 6). By utilizing these 
properties, Zp is rewritten as 



Z P /D l =X 3 + Y 3 + Y- 3 + D 2 X 
+X 2 (D 4 Y + D^ 4 Y- 1 ) 



X{D_ 4 Y 2 + D 4 Y~ 2 ) 



+£> 6 y + D_ 6 y- 1 , (B15) 
D 2 = D 2 /D 1 , 5 ±fc = D ±fc e = F*/-Di(* = 4,6). (B16) 



When the quarks are confined (£ — > 0, £ — > 0), the 0th mod- 
ified Bessel function remains finite (Jo — > 1) but others van- 
ish, I n (y) —s- 0(n > 0,n < 0). In this limit, Zp reproduces 
the standard SC-LQCD expression without the Polyakov loop 
effects, 



Di(ri,fj)= Y, eJW Ai,<, (B8) 

n— — oo 

+ ^ +1 /„_ 2 + ^_ 1 /„ +2 , (B9) 

Ai,2 = — 2(/„ — I n +2lnIn-2) + 5/„+l7„ J„_l 
— 3(7 re+1 J n _2 + In-lIn+2) — In+sInIn-3 
+ In-lIn-2ln+3 + 7n+lA+2 A-3 , (BIO) 

A.± 3 = A T i,i, (Bll) 



7 n=F iA±l + 7„±27 n=F i/ Il z F 2 — 7 Jl ± 2 / n /, lT3 



— 7 n ±l7nJ n =F2 j 
7?n,±5 = A=Fl,=p4 j 

A. ±6 = — (7 n 7„^i + I, 2 l±1 /„ T 3) 

2 

>T 



(B12) 
(B13) 



2(7„ T x^n±l — 7„±27 nT l/ n z F 2 + I n ±2lnln^i) 



(B14) 



Z P -+X 6 - 2X 
sinh((7V c - 



- Y 3 + Y~ 
l)E q /T) 



At 



smh(E q /T) 
•1, A -> -2 . 



2cosh(7V c /i/T) 



(B17) 
(B18) 



In this main part, we redefine Lq = 

5_4, A = A, L 3 = D 6 ,L 3 = 5_ 6 . 



A, A = A, A 
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